Ground state of N=Z doubly closed shell nuclei in CBF theory 



A. Fabrocini^'', F. Arias de Saavedra^-* , CCo'^-', and P-Folgarait^'^-* 
^) Istituto Nazionale di Fisica Nucleare, Dipartimento di Fisica, Universitd di Pisa, 1-56100 Pisa, 

Italy 

^) Departamento de Fisica Moderna, Universidad de Granada, E-18071 Granada, Spain 
. ) Istituto Nazionale di Fisica Nucleare, Dipartimento di Fisica, Universitd di Lecce, 1-73100 

a^ ■ 

Lecce, Italy 

O . 

D . ) Centra Ricerche e Sviluppo, Gruppo Lucchini, 1-57025 Piombino, Italy 

Q 
o 



Abstract 
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The ground state properties of N=Z, doubly closed shell nuclei are studied 

l> 

within Correlated Basis Function theory. A truncated version of the Urbana 
vi4 realistic potential, with spin, isospin and tensor components, is adopted, 

Ij . 

^ . together with state dependent correlations. Fermi Hypernetted Chain inte- 

gral equations are used to evaluate density, distribution function and ground 
^ . state energy of ^^O and ^''Ca. The nuclear matter Single Operator Chain 

approximation is extended to finite nuclear systems, to deal with the non 
commuting part of the correlation operators. The results favourably compare 
with the variational Monte Carlo estimates, when available, and provide a 
first substantial check of the accuracy of the cluster summation method for 
state dependent correlations. We achieve in finite nuclei a treatment of non 
central interactions and correlations having, at least, the same level of ac- 
curacy as in nuclear matter. This opens the way for a microscopic study of 
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medium heavy nuclei ground state using present days realistic hamiltonians. 



Typeset using REVT^iX 



It is now a largely accepted fact that the wave function of strongly interacting nuclear 
systems shows large deviations from independent particle models (IPM). These effects may 
be ascribed to the presence of correlations between the nucleons, coming from the nuclear 
interaction. Several nucleon-nucleon (NN) potentials are presently available, all of them 
fitting the deuteron and the NN scattering data up to energies of several hundred MeV. 
However, their complicated structure and dependence on the state of the interacting nucleons 
has severely hindered the achievement of realistic, ab initio studies of most of the nuclear 
systems. 

The situation is satisfactory for light nuclei. Faddeevi, Green Function Monte Carlo 
(GFMC)i and Correlated Hyperspherical Harmonics Expansion (CHHE)i theories solve 
exactly the Schrodinger equation in the A=3,4 cases for realistic hamiltonians. Recently 
GFMC has been extended up to A=7§. Moreover, these theories (particularly Faddeev and 
CHHE) are now succesfuUy used to study low energy reactions involving three nucleonsi. 
Light nuclei properties may be also described by variational Monte Carlo (VMC)il methods. 
If the spanned variational wave function space is large enough, then the description provided 
by a variational approach is quite accurate (even if not exact). One of the major advantages 
of VMC is its larger flexibility, resulting in the possibility of the extension to heavier nuclei, 



as 



16 



00. 



At the opposite asymptotic side of the nuclear table, infinite nuclear matter has attracted 
the attention of the researchers, as it is thought to be a reliable model for the interior of 
nuclei. High density neutron matter and asymmetric nuclear matter are also objects of 
intensive investigations because of their astrophysical relevance. The equations of state 
(EOS) of infinite systems of nucleons have been studied, in non relativistic approaches 
and using realistic interactions, either by Brueckner Bethe Goldstone (BBC) perturbation 
theoryiil or Correlated Basis Function (CBF) theory0'0. These theories give consistent 
results at densities close to the nuclear matter empirical saturation density {pnm = 0.16 
fm~^), whereas large discrepancies appear at higher density values. A question still to be 
answered is the convergence of the hole lines expansion, on which BBC is based, in the case 



of the continuous choice for the auxihary potentiaL Recent BBG results are all obtained 
within the two-hole line approximationll. Attempts are under way to evaluate the three-hole 
line contribution^. 

CBF nucleonic EOSs give a good microscopic description of nuclear matter around sat- 
uration and provide a description of the neutron stars structure in agreement with current 
observational data0. Moreover, nuclear matter dynamical quantities, as electromagnetic 
responsesEi'0 and one-body Green functionslii, may be accurately addressed by CBF based 
perturbative expansions. 

Medium-heavy nuclei still lack microscopic studies with realistic hamiltonians. In a series 
of papers, the authors succeeded in extending CBF approaches to the ground state of doubly 
closed shell nuclei (both in Is and jj coupling) with semirealistic, central interactions and 
simple two-body correlations, depending only on the interparticle distances and, at most, on 
the isospin of the correlated pairEHll. Nuclei ranging from ^He to ^°^Pb were investigated in 
those papers by model hamiltonians. Aim of the present work is to extend those studies to 
NN interactions and correlations containing spin, isospin and tensor components. We shall 
consider ^^O and "^"Ca nuclei, having doubly closed shells in Is coupling. We shall adapt to 
these systems the cluster summation technique, used in symmetric nuclear matter for state 
dependent correlations. Modern interactions have also important spin-orbit parts, that are 
not included in the present treatment, as well as other remaining components. They will be 
objects of future works. First order cluster expansion has been recently used to study the 
influence of state dependent correlations on one-body density matrix of closed shell nucleiEl. 

Our work is carried out in the framework of the non relativistic description of the atomic 
nucleus with hamiltonians of the type: 

i i<j i<j<k 

The two- and three-nucleon potentials, Vij and Vijk, are determined at large interparticle 
distances by meson exchange processes. The intermediate and short distances regime is 
usually treated in a semi-microscopic or purely phenomenological fashion. We shall use a 



truncated version of the realistic Urbana f 14 model (U14) of NN interactions but we shall 
not consider three nucleon potentials. We shall also present results for the central semi- 
realistic interaction S3 by Afnan and Tang, f^s^, which reproduces the s-wave two-body 
scattering data up to roughly 60 MeV and gives values of the ground state properties of 
light nuclei and of nuclear matter close to those obtained by more realistic interactions. The 
S3 potential has been supplemented in the odd channels, where it is not defined, with the 
repulsive term of the even channels. 

The full U14 has the following parametrization 

^14., = E ^'ir,,)O^P (2) 

p=l,14 

with 

® [1, r, ■ Tj\ , (3) 

Sij = (3f.y ■ (Tjfjj ■ cTj — cTj ■ cTj) being the usual tensor operator. In the fg truncation we 
shall retain components up to the tensor ones, so neglecting the spin-orbit and higher terms 
(p > 6). S3 does not have the p = 3,6 tensor parts. 

The ground state correlated A-body wave function is given, in our CBF approach, by 

^(1,2...A)= (sY[F,,}j^l,2...A), (4) 

where a symmetrized product of two-body correlation operators, Fij, acts on the mean 
field state, $(1, 2. ..A), taken as a shell model wave function built up with (pai'i) single particle 
wave functions. Consistently with the interaction, is chosen of the form: 

F^, = E nr^,)Of,. (5) 

p=l,6 

The tensor components of F^j are omitted in the S3 case. 

The /^(r) functions contain a set of variational parameters determined by minimizing 
the ground state expectation value of the hamiltonian, (H) = (\l'|if|^)/(\l'|\l'). The many- 
body integrals needed for the evaluation of (H) , as well as of the expectation value of other 
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operators, could be in principle sampled by Monte Carlo (MC) techniques. However, MC 
methods for realistic, state dependent models can be efficiently used only in light nuclei. 
An alternative approach, suitable to heavier systems, is the cluster expansion and Fermi 
Hypernetted Chain (FHNC)ii integral equations. FHNC allows for summing infinite classes 
of Mayer-like diagrams and it has been widely applied to both finite and infinite, interacting 
systems with state independent (or Jastrow) correlations . 

The strong state dependence of Fjj, needed for a reahstic description of nuclear systems 
and resulting in the non commutativity of the correlation operators, prevents from the 
development of a complete FHNC theory for the correlated wave function of eq.(|) and 
forces to look for suitable approximations. The single operator chain (SOC) approximation 
presented in ref.(^) (hereafter denoted as PW) for the operatorial correlations, in conjuction 
with a full FHNC treatment of the scalar part, provides an apparently accurate description 
of infinite nuclear and neutron matter^. FHNC/SOC is thought to effectively include the 
contribution of many-body correlated clusters. However, no exact check for FHNC/SOC 
is presently available in nuclear matter, apart the evaluation of some additional classes of 
diagrams. The estimated accuracy in the ground state energy has been set to less than 1 
MeV at saturation densityii0. 

We shall use FHNC/SOC theory to study the ground state of ^^O and '^°Ca described 
by the correlated wave function (^. The ^^O results will be compared with the MC calcu- 
lations of ref.(§), where the scalar part of the correlation is exactly treated by MC, and the 
contribution of the operatorial components (p > 1) is approximated by considering up to 
four-body cluster terms. Higher order contributions are then extrapolated. 

The plan of the paper is the following: in section 2 we present the FHNC/SOC theory 
for the one-body density and the two-body distribution function; the ground state energy 
calculation is discussed in section 3; the results obtained for ^^O and ^''Ca are presented and 
discussed in section 4; conclusions and future perspectives are given in section 5. 
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I. THE FHNC/SOC THEORY FOR FINITE SYSTEMS. 



In discussing the FHNC/SOC approach to the one- and two-body densities, pi(ri) (OBD) 
and p2(i"i)J"2) (TBD), defined as 

Pi(ri) = (E'^(ri-r.)) (6) 

i 

and 

P!(ri,r2) = (E'^(ri " ^i)K^2 - r,)Of^), (7) 

we shall heavily rely on the formalism developed in PW and in ref.(lii), denoted as COl 
hereafter. Most of the quantities we shall introduce and use in this section are described 
in those papers and will not be discussed here. Moreover, the various p = 1,6 components 
of the correlation (and other quantities) will be often referred to as c (p = 1) and, with an 
obvious notation, as a, t and t (tensor). 

In Jastrow FHNC theory, the TBD is written in terms of the scalar correlation, f'^{r), 
and of the nodal (or chain) and elementary (or bridge) functions, N^yivi, Y2) and E^yiri, r2), 
representing the sums of the diagrams having those topological structures, respectively. The 
diagrams are further classified according to the exchange character, (xy), of the external 
points (1,2), xiy) = d,e with d =direct, e =exchange. (cc) (c =cyclic) diagrams are also 
present, whose external points are joined by a single, non closed exchange loop (see COl for 
more details). 

When operatorial correlations are introduced, the nodal functions become N^y{ri,r2), 
where p denotes the state dependence associated with the function. A complete FHNC 
treatment for the full, state dependent TBD is not presently possible, so SOC approximation 
was introduced in PW. It consists in summing p > 1 chains, where each link may contain 
just one operatorial element and central dressings at all orders. We recall that operatorial 
dependence comes also from the exchanges of two nucleons. In fact, to every exchange lines 
forming a closed loop, but one, is associated the exchange operator Pij = — J2q=c,(7,T,aT Olj/'i. 
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The FHNC/SOC integral equations for A^|j^(ri, r2), with x{y) — d,e, are: 
^^.(1. '^) = Y.Y.j d'rsCTs^XUh 3)K%.(3) [x;,^(3, 2) + 7V;,^(3, 2) 

x'y' qr 



(8) 



The allowed {x'y') combinations are: dd, de and eo?, and the coordinate Fj is indicated as 
i. V^J'y,{?)) are vertex corrections in rs that will be discussed later; ^fgf are angular couplings 
given in PW (eqs. 5.6-5.11). Actually, they were given only in the operatorial channels 
(p, g, r > 1). In the p — 1 channel, the coupling function is one if g = r = 1, otherwise it is 
zero. The Xly{l, 2) hnks are defined in COl, while, for p > 1, we have 



Xyi,2) = /i^'(l,2)/i'^(l,2)-7V|,(l,2), 



(9) 



XI{1, 2) = h'^il, 2) {/.^(l, 2)iVJ,(l, 2) + [nn,)f Nlil, 2)} - Nl{l, 2), (10) 



X^l, 2)= h%l, 2) {/^^'(l, 2) (N^l, 2)N:,(1, 2) + 7V,^,(1, 2)) + [^n^)]' (iVfe(l, 2) 
+ Nlil, 2)7V|,(l, 2) + 7VJ^(1, 2)7Vf,(l, 2) - 4 (Ar^^,(l, 2) - po(l, 2)f A^) } 
-iVf,(l,2). 



(11) 



with 



h^{l, 2) = r(ri2) {2r{ru) + nn,)N',,{l, 2)} , 



(12) 



h'^{l, 2) = exp [N^{1, 2)] and = 1 for p = c, u, r, err, otherwise is zero. po(l, 2) is the 
IPM density matrix, given in COl, N^^{1, 2) is the central cc nodal function, given later, 
and Xg^(l,2) = X^^{2,1). The FHNC/0 approximation (corresponding to set to zero the 
elementary diagrams) has been assumed in the above equations. Its validity will be discussed 
in the results section. 

For the cc-type nodals we have 



qr 



KAh 2)= E / ^'^3eSf^c^c(l, 3)Kf (3) fc(3, 2) + KM 2) + N;M, 2) 



(13) 



+ E / d'rs^tE A'^X^Jl, 3)Kf (3) [x^M 2) + A^L(3, 2) + N;^, 2) 

qr>l 
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NU^, 2)= / d'rs^!EX^,,{l, 3)VZ{S) [-po(3, 2) + N^,^(3, 2) + N;^{3, 2) 



+ 



E / rfVaCSf A^X,^,(1, 3)Kf (3) [K^iS, 2) + 7V;,(3, 2) 



qr>l 



(14) 



^^p='(l,2)= / rfV3[-po(l,3)]Kc'(3)A^.y3,2) 



d'r,[-po{l,3)] Kc(3)-1 -po(3,2) + iV; ,(3,2) 



(15) 



A^^P>'(1, 2)= E / ^'^3e?3^ [-Po(l, 3)A^] Kr(3)A^.^p(3, 2) (16) 
+ E / ciVaCfSf [-po(l,3)A''] [Kf (3) - l]iV;,(3,2). 

qr •' 

Again, ^^^(1,2) is defined in COl, and 

XI>\1, 2) = hF{l, 2)h'^{l, 2) [iV,^,(l, 2) - po(l, 2)] + {[r(ri2)]' /i^(l, 2) - l} N^,,{1, 2). (17) 

The x(p) subscript indicates that the external point is (not) reached by a X hnk, N^^ — 
+ + + N^^ and TV^Jl, 2) = Arp^(2, 1). 

Because in an exchange loop involving more than two nucleons only one of the exchanged 
pairs does not have any operatorial link from and in the spirit of the SOC approximation, 
the correlations appear once in the cc chains, just for that pair. 

Single operator rings (SOR) were also approximately included into the central chains in 
PW. SOR are closed loops of operators having a non zero C-part. A product of operators 
can be expressed, by the Pauli identity, as the sum of a spin and isospin independent piece 
(the C-part) and a remainder, linear in or Tj. In nuclear matter, as well as in the nuclei 
we are considering (doubly closed shell nuclei in Is coupling), the spin- isospin trace of the 
remainder vanishes leaving only the C-part as the final contribution of the product. An 
example is 0^2^23^31! having a C-part of 3. SOR contribute to p = 1 chains and were 
introduced by PW in the definition of X^^. A drawback of this approach is that touching 
SOR, i.e. SOR having a common vertex, are wrongly counted because commutators are 
neglected. For this reason, we do not follow PW and do not include SOR in our treatment 
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of N^y. However, in one test case we have used the PW prescription to gauge the relevance 
of the missing contribution. The results will be presented later. 

The OBD, pi(ri), is computed following COl. Its structure results to be 

Pi(ri) = P^(ri) [1 + C/7(ri)] + U:^{r,)C,{r,), (18) 

with 

p^(rO = [po(ri) + ^e^(rO]Q(rO, (19) 

where po(ri) = Ea|0a(l)P is the IPM density, Crf(ri) = exp {^/^^(ri)} and t/rf%)(ri) 
are the central vertex corrections of COl. They are solutions of two integral equations, 
(A. 19) and (A. 20), given in the appendix of that paper. These equations must be modified 
because of the presence of operatorial correlations by the substitutions ^d{T^2) Pi(r2) and 

The SOC operatorial vertex corrections, C^^('e)(ri), are solutions of the equations 

E / ^'^2/^^(1, 2)r{r^,) {[Fin,) + nrn)N^^{l, 2)] pl{2) (20) 
+ [nn2) + 2r(ri2) Ar^,(l, 2)] N^l. 2)C,(2) 
+ r(n2)A^fe(l,2)Q(2)}, 

C/f (1)= I d'r,h\l, 2)/^(ri2) {[/^(ris) + 2/'=(ri2)Ar,^,(l, 2)] (21) 

p>i 

[Nl,{l, 2)p1{2) + {N^,,{1, 2)N:,{1, 2) + Ar|,(l, 2)) C,{2)] 

+ mru) [Keih 2)Q(2) + AT^^l, 2) (p^(2) + Arj,(l, 2)C,(2))]} 

where Ap='^'^ = 1, 3, 6, 3, 9, 18 and 

^7(1) = -8 E A"^' I / d'r2h^i^, 2)/^(ri2)r(ri2) [-po(l, 2) + Kc{l, 2)f C,{2) (22) 
+ j cf'r, J dV3^,^,(l, 2)^,^,(1, 3)/i^(2,3)/nr23)r(r23) [-po(2, 3) + Ar^^,(2, 3)] C,(2)C,(3)| , 
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with g'^^Xi.j) = [nr,,)fh%i,j) [-p,(^,J)+N^J^,J)]. 

The vertex corrections of the nodal equations, V^y, can be expressed in terms of the 
OBD and of the [/-functions. For the central, p = 1 chains we have V^^~^^{i) = pi{i) and 
^de7d\ecJS) — Cd{i) [1 + U^^{i)]. As far as the p > 1 chains are concerned, we insert central 
vertex corrections only : V^^{i) = pl{i) and VX^ed.ee.cclO = Cd{i)- 

II. ENERGY EXPECTATION VALUE. 

In order to evaluate {H), we use the Jackson- Feenberg identity^ for the kinetic energy, 
with the result 

(T) = Tjp = T^ + Tf, (23) 

with 

Ts = -A — ($*G'V?$ - (Vi$*) (Vi$)) (24) 
4m 

and 

Tf = -A — ($* GV?G-ViG-ViG' $), (25) 
where G = SY[Fij. In turn, is written as 

T^ = T^'^+Ti'^+T!^'K (26) 

The T^p""^ terms correspond to contributions where the kinetic energy operator acts on a 
nucleon not involved in any exchange {n = 1) or belonging to a two-body {n = 2) or to a 
many-body {n = 3) exchange loop. 

For T^^ we obtain 

Ti'^ = -^J d'r,pr,{r,)CM) [1 + ^^7(^1)] (27) 
and pri(ri) is given in COl. 
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For the remaining parts of Tjp, as well as for the two-body potential energy (v) — V2, 
we are faced with computing the expectation values of specific two-body operators (apart a 
small three-body term, in which will be discussed separately). 

We start with Tp + V2 — W, also called interaction energy in PW, and define //j^(ri2) 

as 

Hj%u) = {r(ri2)VV'(ri2) - Vfiru) ■ Vfir,^)} + f{ruy{ru)f{ru). (28) 

In FHNC/SOC, W is split into four parts 

W^Wo + Ws + Wc + Wcs, (29) 

where Wq is the sum of the diagrams with only central chains between the interacting 
points (IP), connected by Hjp. Ws sums diagrams having SOR touching the IPs and central 
chains; Wc contains diagrams with one SOC between the IP and Wcs contains both SOR at 
the IP and SOC between them. 

Wq is given by 

Wo= I J d'r, J d'r2Hy^{r,2W{l, 2) {k^^'A' [p^(1)p^(2) + pl{l)C,{2)NUl, 2) (30) 
+ C,il)pl{2)N^,,{l, 2) + C,(1)C,(2) {N:,{1, 2) + 7V^^,(1, 2)7V^^,(1, 2))] 
_ 4ii'^^-^i^*™A™A™Cd(l)Cd(2) {N'Jl, 2) - po(l, 2))'} . 

A sum over repeated indeces is understood and the matrix K^^'^ is given in PW. 

The presence of Ws is due to the non commutativity of the correlations. In nuclear 
matter and for state independent correlations, this term is absent because of the complete 
cancellation of the separable diagrams (see PW for a more complete discussion of this point). 
We obtain 

Ws^ \ I d'r, J d'r^ny^irumi, 2) [k^^'A' (l + Re±2£±^^ (3I) 
'{pl{2) + C,(2)iV^,(l, 2)) (p?(l)M^(l) + C,(1)M^(1)) 
+ imKaih 2) + C,{2) (iV,^,(l, 2) + iV,^,(l, 2)iVJ,(l, 2))) C,(1)M^(1) 
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+ (p^(2) + C,(2)iV,^,(l,2))C,(l) (2M^i(l) + M^2(i; 



Din + D 



1 + 



+ 1^2} 



2a 



In 
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C,(1)C,(2) [iV,^,(l,2)-po(l,2)]'M,"(l) 



are given in eq.(5.23) of PW and we consider only terms linear in the vertex 
corrections, taken of the simplified form 



1 = I dh, firu) ■ h'^il, 2) {pli2) + C,(2)iV^,(l, 2)} , 



(32) 



MJ>'= A'jd'r^ [f{n,)\ h%l,2) {pI{2)N:,{1,2) 
+ C,i2) [iV:,(l, 2) + iV,^,(l, 2)iV:,(l, 2)]} , 



(33) 



/>i 

cl 



-4A'A' I d\2f{n2)f{ruW{l, 2)C,{2) [iV:,(l, 2) - po(l, 2)]' , (34) 



= -4A' / /'(ri2) "/i'=(l,2)C,(2) [iV,^,(l,2) -po(l,2)]^ 



(35) 



Contributions from SOCs have not been inserted into M^. In Wc we must keep track of 
the order of the operators both in the IP and in the SOC. Its expression is quite lenghty and 
it is given in the appendix. Because of the large number of involved operatos, Wcs is the 
messiest term among all, and the smallest, so we approximate W^^ ~ /^^o]) where 

the j-index refers to the j-component of i^xF- A more involved factorization approximation 
to Wcs was used in PW and its validity was set to within ~ 0.2 MeV. We have checked our 
approximation against the PW one in nuclear matter and found agreement up to the second 
decimal digit. 

The decomposition p9|) can be carried on also for TP-* and Tjf'"^^ (the two-body part of 



tP'*). The result is: 



T% = d'n I f/V2PT2(l,2)C,(l) {A-'=U'A'Crf(2) [rin2)fHri2Wil,2 



(36) 



5^l5kl] + Crf(2) - 1} 
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(3,2) 



0,0 



-2- 



m 



d'r, j dV2PT3(l,2)C,(l){ir*'='AWQ(2)iV,^,(l,2) 



f(n2)/'(ri2)/^^(l,2)-5,i4i +Q(2) iV,^,,(l,2) + iV;,(l,2) 



+ (C,(2)-1) [iV,yi,2) + iV;^(l,2)]} 



(37) 



Tf} = --j d'nj £r,pT2il,2)Ca{l) i^K^'^A'A' (l + 

C,{2) [r (ri2)f (ri2)/^^(l, 2) - 5,,5ui] {M^{1) + M7(2)) 
+ t/7(l) (C.(2)-l) + C.(l)f/7(2)}, 



fcm 



(38) 



T. 



(3,2) 



d'n / c?Spt3(1, 2)C,(1) ir^'^'A'A' 1 + 



Dkm + D, 



Im 



(39) 



C,(2)iv:,(i, 2) [f{n,)f{n,)h%l, 2) - (M™(1) + M™(2)) 

+ c,(2)iv^^,(i, 2) ([/7(i) + [77(2)) - f/7(i) [n:^{i, 2) + iv;^(i, 2)] } . 

(2) (3 2) 

pT2,z are defined in COl. Again we give T^ ,: and T^_^ in the Appendix and the cs terms 
are evaluated according to the approximation used for W^s- 

A three-body term, t'^''^\ originates from Vi$* ■ Vi$ in eg . (^4|) . It was not computed 
in COl because it is known to provide a small contribution in nuclear matteiS. Here we 
include it, adopting the Kirkwood superposition approximation (KSA) for the three-body 
distribution functionS and considering only the p = 1 correlations contribution. Following 
these prescriptions, we obtain 



T?-'^ = -2'— j d'n j d'r2 j d'rsVipoil, 2) ■ Vipo(l, 3)C,(1) (40) 
{C,{2)C,{3)g^U^, 2)g^,,{l, 3) M,{2, 3) - 1) {K,{2, 3) - po(2, 3)) + Ar,^,(2, 3)] 
+ [C,(2)(?^,(l, 2)Ar,^^(2, 3) (^7^,(1, 3)C,(3) - 1) + 2 ^ 3 
+ [Ar;^(2,3) - po(2,3)] [^7^,(1,3)^(3) - 1] [^7^,(1, 2)C,(2) - 1]} , 



where gddihJ) = [f%r,j)]^h%i,j) . 



14 



III. RESULTS. 



All the results presented in this section have been obtained with the single particle wave 



functions, 4>a{i), generated by a harmonic oscillator well with oscillator length b = \Jfi/mui. 
In principle, h could be considered as a variational parameter; however we kept it fixed, at 
h = 1.543 fm for ^^O and b = 1.654 fm for ^°Ca, because our aim here is to develop and assess 
the finite nuclei FHNC/SOC theory, rather than to perform a fully variational calculation, 
to be compared with experimental data. This problem will be takled when the complete, 
realistic hamiltonian will be within reach of our approach. 

The best choice for the correlation operator Fij is obtained by the free minimization of 
the FHNC/SOC energy functional and the solution of the corresponding Euler equations, 
[S{H) /6Fij = 0]. This method is not practicable for realistic NN potentials and one has to 
resort to less ambitious ones. In COl two types of correlations were investigated: a simple, 
two parameters gaussian, fcif'), and a more effective i?n/er function, fEui{f), obtained by 
minimizing the energy evaluated at the lowest order of the cluster expansion, {H2). The 
latter was also adopted in the nuclear matter studies of PW. Here we shall use the Euler 
correlation corresponding to eg.®, so extending to the state dependent, finite system case 
the approach of COl and PW. 

Without going into many details, the correlation is computed in the (T, a) channels, with 
a = {S,t), where T and S denote the total isospin and spin of the pair and t the tensor 
part (5=1 for the t-channel). In the 5* = case, /to(^) is solution of the Schrodinger-like 
equation 




V^FTo(ri2) + Vroiru) - Xto i^To('"i2) = 0, 



(41) 



m 



where Frairu) = f t a{r 12) Prsi.^ 12)^ 



Pt5(1, 2) = Po(l)po(2) - 16pg(l, 2)(-)^+^ 
1 

Qt5(1,2) = -VTs{ri2)PTs{l,'2) + _p^2(l,2)(-) 



T+S 



(43) 



(42) 



and 
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2^ 



^T5(n2), with X = {Q,P), is defined as in eq.(3.10) of COl. 
The S = 1 correlations are solutions of two coupled equations 

V^Friiru) + Vniru) - \ti Friir^) + 8 [vrtiru] - Xxt] Frtin^) = 0, (45) 



m 



and 



2 



V'Frf(ri2) + 



6 

VTi(n2) - 2t;Tt(ri2) H + ^^t* - Ati 



FTt(ri2) (46) 



+ [vTt{ri2) - Xrt] Friiru) = 0. 



These equations are solved under the healing conditions frsi^ > (Its) = 1, frtir > 
drt) = and /■^Q,(r = dxa) = 0, where dxa play the role of variational parameters. The 
/^(r) correlation functions are then obtained by fxaij) (see PW). 

In nuclear matter only two healing distances are used: dc and dt-, for the central {dTs) 
and tensor (rfrt) channels, respectively. We make here the same choice. Additional nuclear 
matter variational parameters are the quenching factors of the NN potentials in the Euler 
equations. As in PW, we take = 1 and a^^^ = a. We have already stated that, for the 
time being, it is not our interest a full variational search, so we have taken the nuclear matter 
parameters given in ref.(0) for U14. They are: dc = 2.15 fm, dt = 3.43 fm and a = 0.8. 

The '^'^Ca correlation functions are shown in Figure 1 and compared with the correspond- 
ing nuclear matter functions, at saturation density. They are similar, especially the longer 
ranged tensor ones. The most visible differences are found in the a and r components and 
in the shortest range part of f^. We stress that additional differences could arise from the 
minimization process, as the energy mimimum will probably correspond to a different choice 
of the parameters. 

A measure of the accuracy of the FHNC/SOC approximation is how well the densities 
normalization sum rules are satisfied: 
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Si = J d^np,{r,) = A, (47) 
^2 = J^X^) I I ^'^2P2(ri, r^) = 1, (48) 

Sr = ^j d'n I d\2Pl{vi, V2) = -1. (49) 

The spin saturation sum rule, = —1, holds only in absence of tensor correlations^. 
Both the TBD and its sum rules are evaluated following the decomposition (^1^. 

Deviations of the sum rules from their exact values are due to {i) the FHNC/0 scheme 
and (ii) the SOC approximation. The influence of the elementary diagrams was addressed 
in COl. It was found that E^^'^^, i.e. the sum of the ee-elementary diagrams whose external 
points belong to the same exchange loop, may substantially contribute to both Sr and to 
the potential energy, if the potential has large exchange terms. This fact can be understood 
if we consider that a four-point elementary diagram, linear in the central link, [f'^]'^ — 1, is 
contained in E^^'^^, as well as diagrams linear in the operatorial link f^f^^^. The insertion 
of these diagrams in the FHNC equations was termed as FHNC-1 approximation, and we 
shall keep this terminology. 

Results for the sum rules are presented in Table I for different models of the correlation: 

(Jastrow , p = 1 component only), and (without and with tensor correlations, 
respectively). The Table shows also the FHNC-1 corrections. In all cases, 5*1 shows a largest 
error of less than 1%. This is also the accuracy that we find in the Jastrow case for S'2,t, in 
FHNC-1, as already noticed in COl. The situation is worst for the operatorial correlations, 
where FHNC/SOC violates the sum rules by a maximum amount of ~ 9%, similar to what 
was already found in nuclear matter in ref.(0). 

The ground state energetics is displayed in Tables II-IV for ^^O and Tables V-VII for 
^°Ca, for each correlation model. The columns (0, s, c, cs) show the contributions to the 
two-body operators expectation values, as given in eg. (^9]) . The ^^O results are compared 
with the calculations of PiepeiS with the cluster Monte Carlo (CMC) method of ref. (i). As 
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already outlined, in the CMC approach the Jastrow part of the correlation is exactly treated 
by MC sampling and the remaining operatorial contributions are evaluated by MC up to 
four nucleon clusters. The fourth order cluster expansion seems to be enough to provide 
a rehable convergence and to consider the CMC numbers as a benchmark for FHNC. The 
Tables also contain the FHNC-1 corrections. 

The ground state energy mean value, Egg, is then given by Egg = {H) — Tcm, where Tcm 
is the center of mass kinetic energy, whose calculation is discussed in COl. 

For the Jastrow correlation, FHNC shows an error of 1 — 2% for (T) and {V2) in ^^O. 
The total energy percentile error is bigger (~ 9%) as {H) is given by the cancellation of 
two large numbers. We meet the same situation in the and models. The kinetic and 
potential energy errors are 3 — 4% in the first case and 5 — 7% in the latter. The absolute 
error in (H) is well less than 1 MeV in all models. Again, this finding is consistent with the 
estimated accuracy of FHNC/SOC in nuclear matter. We notice that most of the binding 
is given by the OPE parts of the potential, v""^ and v*"^. In absence of the last component, 
^^O is not bound in our model. The same holds in ^'^Ca, where the introduction of tensor 
correlations and potentials increases the kinetic energy by ~ 5.6 MeV, compensated by an 
additional potential energy contribution of ~ —13.6 MeV, providing a bound nucleus in 
the case. For the sake of curiosity, we recall that the experimental binding energies per 
nucleon are —7.72 MeV in ^^O and —8.30 MeV in ^°Ca, to be compared with the computed 
values Egg = -5.15 MeV (^^O) and -7.87 MeV (^°Ca). 

In Table Vlll we compare the expectation values of the components of the potential 
with the nuclear matter results, within the same FHNC/SOC approximation and in the 
model. It is interesting to notice a kind of convergence with A for the potential energies, in 
particular for the large OPE related components, whose contributions, in ^°Ca, are already 
very close to the nuclear matter values. We stress that a more meaningful comparison would 
imply the use of Hartree-Fock single particle wave functions, or, at least, a minimization on 
the single particle potential parameters. 

In Table IX we show the influence of the SOR in ^^O. SOR have been inserted according 
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to PW. In general, they contribute for less than 1% of the FHNC/SOC value, with the 
exception of (f'^), where they give a 17 — 18% contribution, actually worsening the agreement 



The effects of the correlations on the ground state structure are shown in Figures 2 and 
3, giving the OBDs and the two particle distribution functions, P2(i"i2), defined as 



where R12 = | (ri + is the center of mass coordinate. In both figures, the the 
Jastrow and the independent particle models are compared 

Large parte of the reduction respect to the IPM is due to the Jastrow, short range 
correlations. The operatorial correlations slightly enhance the OBDs, as in the first order 
cluster analysis of ref.(lli). They have the same effect in P2(ri2)5 where the dip at short 
distances is due to the repulsive core in the nuclear interaction, as already found for A=3,4 
nuclei in ref.(c2l). 




At the beginning of this section we have explained why we did not look for a variational 
minimum for the truncated version of U14. However, it is certainly of interest to try to 
understand how reliable are the nuclear matter parameter values and how far they are from 
the true minimum. To this aim we have minimized, with respect to dc-, the energy for the 
S3 model described in the Introduction (keeping the same harmonic oscillator wells as U14). 
The results are displayed in Table X. The first row corresponds to the U14 nuclear matter 
dc = 2.15 fm, whereas the second gives the computed minima. The minimization produces 
a small gain in the binding energy and S3 appears to provide two nuclei underbound of ~ 1 



In this article the FHNC technology developed in COl to describe finite nuclear sistems 
has been extended to state dependent correlations containing up to tensor components. As 



with CMC. 




(50) 




MeV. 



IV. CONCLUSIONS. 
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in infinite nucleon matter, the non commutativity of the two body correlation operators 
does not allow for a complete FHNC treatement, which is instead possible for purely scalar, 
Jastrow-type correlations. The single operator chain approximation scheme (which is ef- 
fectively employed in nuclear and neutron matter) has been extended to the finite case. 
The resulting set of integral equations has been solved either by neglecting the class of 
the elementary diagrams (FHNC/0 approximation) or by considering only the lowest or- 
der elementary contribution in the dynamical correlation (FHNC-1). As an application, we 
have studied the ground state properties of the doubly closed shell nuclei in the Is coupling 
scheme, ^^O and "^^Ca, interacting by the central and tensor components of the realistic 
Urbana vu nucleon-nucleon potential. 

The analysis of the sum rules shows that the FHNC/SOC equations provide a consid- 
erably accurate one-body density, whose normalization is violated by much less than 1 %. 
A comparably good accuracy is obtained for the normalization of the central component of 
the two-body density (^'2 = 1 sum rule) when tensor correlations are not included. If they 
are considered, then the excellent fuUfiment of S2 obtained with purely central correlations 
slightly worsens. The inclusion of the first order Jastrow elementary diagram in FHNC-1 
does not improve the outcome. The inclusion of the analogous diagrams, linear in the oper- 
atorial correlations, is presently under consideration. In any case, the worst violation of the 
sum rule is ~ 9 %, close to what was found in nuclear matter. A similar situation is met for 
the isospin saturation, S-r — —1, sum rule. 

The various energy contributions in ^^O have been compared with those obtained within 
the cluster Monte Carlo approach. The maximum disagreement with the CMC results varies 
from ~ 2 % for the Jastrow model to ~ 7 % for the tensor model. The absolute error in the 
ground state energy per nucleon is always well less 1 MeV, compatible with the estimated 
accuracy of the FHNC/SOC approach in nuclear matter at saturation density. 

The same truncated vu interaction has been also used to study the ground state of ^°Ca. 
We have verified that for both ^^O and "^^Ca only the insertion of the long range one pion 
exchange parts of the potential (and related correlations) binds the nuclei. 
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No minimization along the correlation and single particle potential variational parameters 
has been carried on, but we have rather taken the nuclear matter values. We have postponed 
this task to future works, when a completely realistic hamiltonian will be within reach of 
our method. However, a partial minimization on the correlation healing distance, dc, for the 
simpler, central Afnan and Tang potential seems to point to little variation of the parameters 
in going from the infinite to the finite case. 

Even if this is still an intermediate step towards a full microscopic description of inter- 
mediate and heavy nuclei, our results are very promising. In fact, we may conclude that the 
FHNC/SOC approach to finite nuclei shows at least the same degree of accuracy estimated 
in the best variational nuclear matter studies. In this respect, we consider as mandatory the 
inclusion of spin-orbit terms in both the interaction and correlation, as well as the extension 
to the jj coupling scheme, in order to cover all the range of the doubly closed shell nuclei. 
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APPENDIX: 

In this Appendix we give the Wc, T^^j, and T^J, expressions. is given by the sum 

Wc - Wc{dd) + Wc{de) + W^^ed) + Wc{ee) + Wc{cc) (Al) 

with 

WMd)= \ I d'n I dh,W]^{n,)h\l, 2)Ar^>^(l, 2) {[p^(l)p?(2) + p^(l)C.(2)iV,^,(l, 2) (A2) 
+ C,(lK(2)iV,^,(l, 2) + Ca{l)Ca{2) {NI,{1, 2) + Ar,=,(l, 2)Ar,^,(l, 2))] 

-4Q(1)Q(2) [Ar,=,(l, 2) - po(l, 2)]' A" 
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8 ^ ' 
I ^^j^nlm j^ijm' j^mm' k j^k _j_ j^jkmj^mnm' £ilm' _j_ j^nim j^mjm' j^klm'^ ^ 

iye(cie)= lye(erf) = ^ / rfVi I d^^HyXr^mh 2)iVi>^(l, 2) M(1)C,(2) (A3) 



WJee) 



W,{cc)^ ^(-8) I dVi I dh2H%{T,2W{l, 2)C,(1)C,(2) [Ar^^,(l, 2) - po(l, 2)] (A5) 
A'^ {^^SLll, 2)K^^^K^^^'L^'^' + [Ar^>],(l, 2) + <>i(l, 2); 

— ^ J^j km j^mnm' j^im' I j^l i j^jkmj^mnm' jj.m'l i j^jkm j^inm' j^mrn' I 

8 V 

I j^jkmj^inm'jynm'l _j_ Q^j^jkmj^imm'j^nm'l _j_ j^ijm j^nmm' j-^m' kl _j_ j^ijm -^nkm' jyn'ml^^ 

In the last equation, N^^^j^) are cc-nodal functions having the hnk reaching the left 
(right) external point. N^^^J^^ has as an internal link, 
is given by 



A^S(l,2) = A^^>;^(l,2) + A/^>i(l,2), 



ri>i 



l>i 



(A6) 



where A^^^^^^ and A^^^^ are solutions of 



2) = E / d'rsaiiXi^Lih 3)Kf (3) [^c^c(3, 2) + N^^^S, 2) + Ar;,(3, 2)] A^ (A7) 
2) = E / ^'r3e?3Uc^c,L(l, 3)Kr(3) [-po(3, 2) + 7V,^^(3, 2) + 7V;^(3, 2)] A^ (A8) 



and 



^cc>i(l, 2) = h\l, 2)h%l, 2) [Ar^,(l, 2) - po(l, 2)] + {[r(r:2)]' h^{l, 2) - l} AT^,, Jl, 2). (A9) 
For the other functions, we have iV^>]^(l,2) = Ar^>l(2, 1) and N'>1, = Ni>^ - N'>i - 

^^cc,R- 
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Finally, TPj and Tff^ are given by 



Tfl = -- / <fn [ d\2PT2il, 2)C,{l)Cai2)fin2)firi2mi, 2)Nj>\l, 2) A" (AlO) 

777- 



nim rmlk 



+ K 



ikm T nlm 



knm T mli 



12 ' 



{2f (ri2)/'(ri2)/i'^(l, 2)7Vi^^(l, 2)iV,^,(l, 2) 



(All) 



ikm j-jilm _j_ j^knm j^mil j^l _j_ j^knm jynli\ 



— r J<[™''^ J^''^^^ -\- 

8 V 

I / ^j^nlm j^mik ^k , j^knm £j,lm , ^''^'"^ 

19 V 



+ [r (ri2)/'(ri2)/i^(l, 2) - SaSki] [(iVg^l, 2) + A^L':]i(l, 2)) 
4 ^ 



nml 



nim T mkl 



nfem r mii 



1,2)K 



} 



-2- / dV, / (iV2PT3(l,2)Cj(l)^'A'{Q(2) 

+ [C.(2)-l]iV^>L(l,2)}. 



Ar^>J„,(l,2) + Ar;>;„,(l,2) 



i>i 



The L*-''^ matrix is given in PW. 
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TABLES 

TABLE I. ^^O and ^^Ca sum rules for U14 with different correlations. The /"^ line corresponds 
to the Jastrow model; the /^(^) line gives results with (without) tensor correlations. Numbers in 
parentheses are obtained in the FHNC-1 approximation. 



f 16.00 0.998 (1.002) -1.057 (-1.001) 

f 16.03 0.988 (1.001) -0.980 (-0.965) 

f 16.01 1.051 (1.054) -0.943 (-0.930) 

f 40.00 0.999 (1.001) -1.067 (-1.002) 

^^Cs. 40.03 1.005 (1.007) -1.074 (-1.056) 

f 39.86 1.089 (1.091) -0.997 (-0.981) 
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TABLE II. Contributions to the energy per nucleon, in MeV, for ^^O with the truncated U14 
potential, the p = 1, Jastrow part of the Euler correlation and the harmonic oscillator single 
particle wave functions discussed in the text. The E^^'^^ column gives the FHNC-1 correction. 



+s +c +CS +E^^''^ CMC 





14.32 








/OA 




3.90 






/Q OA 




0.73 






7.(3,3) 


0.07 








Tf 




5.59 






(T) 




24.61 




24.33(21) 






0.84 


0.88 


0.93(28) 






1.28 


1.25 


1.27(08) 






2.46 


2.40 


2.43(12) 






-27.34 


-26.59 


-26.24(26) 


{V2) 




-22.76 


-22.07 


-21.56(25) 


{H)/A 




1.78 


2.54 


2.77(09) 




0.82 








E,,./A 






1.72 
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TABLE III. As in Table II, for the correlation model 










+s 


+c 


+CS 


1 jjjexch 


CMC 




14.41 


















4.98 


4.97 


4.97 


4.97 










1.44 


1.44 


0.05 


0.05 






^(3,3) 


0.07 














Tf 




8.09 


8.28 


7.80 


7.79 






(T) 




29.00 


29.17 


27.30 


27.29 




26.15(31) 






3.04 


3.13 


2.43 


2.41 


2.41 


2.72(37) 






2.27 


2.26 


2.07 


2.07 


2.02 


2.07(10) 






2.38 


2.35 


2.35 


2.35 


2.28 


2.40(12) 






-34.42 


-34.56 


-32.82 


-32.81 


-32.10 


-31.76(33) 


{V2) 




-26.73 


-26.82 


-25.96 


-25.98 


-25.39 


-24.58(29) 


{H)/A 




2.26 


2.36 


1.34 


1.31 


1.90 


1.57(09) 



T^/A 0.82 

Eg,/A 1.08 
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TABLE IV. As in Table II, for the f correlation model 










+s 


+c 


+CS 


1 jjjexch 


CMC 




14.75 


















5.04 


4.95 


4.93 


4.93 










1.33 


1.29 


-0.04 


0.00 






^(3,3) 


0.07 














Tf 




11.45 


12.22 


11.46 


11.41 






(T) 




32.63 


33.28 


31.16 


31.16 




29.45(33) 






3.03 


3.31 


2.41 


2.33 


2.33 


2.35(43) 






2.17 


2.13 


2.02 


2.02 


1.97 


2.00(13) 






2.34 


2.26 


2.36 


2.36 


2.29 


2.23(14) 






-33.79 


-34.25 


-32.71 


-32.69 


-32.03 


-30.12(42) 






0.31 


0.30 


0.26 


0.26 


0.25 


0.27(01) 






-11.42 


-11.82 


-10.41 


-10.36 


-10.28 


-9.77(09) 


{V2) 




-37.37 


-38.07 


-36.07 


-36.08 


-35.47 


-33.03(31) 


{H)/A 




-4.74 


-4.80 


-4.91 


-4.92 


-4.33 


-4.59(10) 



T^m/A 0.82 
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TABLE V. As in Table II, for ^°Ca. 



+s +c +CS +E':^''^ 





14.81 






m{2) 




5.00 








1.74 




T-.(3,3) 


0.72 






Tf 




8.29 




(T) 




30.55 








-1.45 


-1.41 






1.60 


1.57 






3.06 


2.99 






-33.30 


-32.40 


{V2) 




-30.09 


-29.26 


{H)/A 




0.47 


1.30 




0.28 






Egs/A 






1.01 



30 





TABLE VI. As 


in Table V, for the 


correlation model 









+s 


+c 


+CS 


^E^xch 




14.92 












6.06 


6.06 


6.06 


6.06 






2.74 


2.75 


1.05 


1.04 




^(3,3) 


0.73 










Tf 


11.72 


12.08 


11.39 


11.37 




(T) 


36.17 


36.54 


34.15 


34.12 






1.16 


1.21 


0.18 


0.14 


0.14 




2.78 


2.77 


2.45 


2.45 


2.38 




2.94 


2.90 


2.83 


2.83 


2.74 




-42.69 


-42.92 


-39.78 


-39.76 


-38.92 


{V2) 


-35.80 


-36.04 


-34.33 


-34.34 


-33.66 


{H)/A 


0.36 


0.49 


-0.18 


-0.22 


0.46 




0.28 










Egs/A 










0.18 
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TABLE VII. As in 


Table V, for the f 


correlation model 









+s 


+c 


+CS 


^E^xch 




15.37 












6.12 


6.07 


6.07 


6.07 






2.50 


2.49 


0.92 


0.93 




^(3,3) 


0.70 










Tf 


16.27 


15.63 


16.66 


16.62 




(T) 


40.96 


42.47 


39.72 


39.69 






1.13 


1.29 


-0.03 


-0.22 


-0.21 




2.61 


2.54 


2.35 


2.35 


2.30 




2.83 


2.69 


2.79 


2.79 


2.71 




-41.41 


-42.27 


-39.54 


-39.48 


-38.73 




0.37 


0.35 


0.30 


0.30 


0.29 




-15.43 


-16.23 


-13.72 


-14.22 


-14.14 


{V2) 


-49.91 


-51.61 


-47.85 


-48.48 


-47.28 


{H)/A 


-8.95 


-9.15 


-8.13 


-8.79 


-7.59 




0.28 




















-7.87 
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TABLE VIII. Breakup of the FHNC/SOC potential energies in MeV per nucleon for ^^O, ^°Ca 



and nuclear matter in the 


model. 








160 


40Ca 


nm 




2.33 


-0.21 


-3.04 




1.97 


2.30 


2.46 




2.29 


2.71 


2.84 




-32.03 


-38.73 


-37.69 




0.25 


0.29 


0.32 




-10.28 


-14.14 


-14.06 


Egs/A 


-5.15 


-7.87 


-13.16 
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TABLE IX. Contributions in MeV to the ^^O energy per nucleon with and without SOR 
insertions and in CMC. 







SOC 


SOC+SOR 




CMC 


J 












\J- } 






07 oc; 










2.41 


2.00 




2 72('37) 


iv") 




2.02 


2.02 




2 07(10) 


\^ / 




2.28 


2.30 




2.40(12) 






-32.10 


-32.01 




-31 76f33l 


J 












IT\ 




^1 Ifi 

O -L . -LU 


O -L .UO 






\^ 1 




2.33 


1.92 




2 35 ('43') 


\" 1 




1.97 


1.96 




2.00(13) 






2.29 


2.30 




2.23(14) 


\^ 1 




-32.03 


-31 81 




-30 1 2('42'l 






0.25 


0.26 




0.27(01) 






-10.28 


-10.21 




-9 77(091 


TABLE X. 


Energies per nucleon (in 


MeV) and central healing 


distances (in 


fm) for the S3 


interaction in the 


model. 








160 


dc 


Egs/A 


40Ca 


dc 


Egs/A 




2.15 


-6.08 




2.15 


-7.37 




1.96 


-6.27 




2.02 


-7.50 
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FIGURES 

FIG. 1. Euler correlation functions in ^°Ca and nuclear matter (nm) at saturation density. In 
the right panel, lines are the '^''Ca correlations and symbols are the nm ones. 

FIG. 2. One body densities in ^^O (left) and ^'^Ca (right). The solid lines correspond to the 
model, the dashed to the Jastrow model and the dot-dashed to the IPM. 

FIG. 3. Two particle distribution functions in ^^O (left) and ^'^Ca (right). As in Figure 2. 
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